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Abstract 

The atmospheric neutrino fluxes observed by underground detectors are not symmetric for rota- 
tions around the vertical axis. The flux is largest (smallest) for the u's traveling toward east (west). 
This asymmetry is the result of the effects of the geomagnetic field on the primary cosmic rays and 
on their showers. The size of the asymmetry is to a good approximation independent from the 
existence of neutrino oscillations, and can in principle be predicted without any ambiguity due to 
existence of unknown new physics. We show that there is an interesting hint of discrepancy between 
the experimental measurement of the east-west asymmetry by Super-Kamiokande and the existing 
theoretical predictions. We argue that the discrepancy is a real effect caused by the neglect, in the 
existing calculations, of the bending of the charged secondary particles (in particular of the muons) 
in the cosmic ray showers. The inclusion of this effect in the calculation gives results in quantitative 
agreeement with the data. We comment on the implications of this result for the study of neutrino 
oscillations. 



1 Introduction 



The precise calculation of the fluxes of atmospheric neutrinos has become a much more 
interesting task since the measurements of Super-Kamiokande ^ and other detectors 
0, ^, ^ ^ have given convincing evidence that u flavor transitions exist and are modify- 
ing the intensity and distorting the angular distributions of the u fluxes. As the field of 
atmospheric neutrino is maturing from the "discovery era" to the era of "precision mea- 
surements" the need for detailed, accurate predictions for the the atmospheric neutrino 
fluxes has clearly grown in importance, since comparison of data and predictions is now 
used to extract from the data the parameters that describe exciting new physics. At the 
same time the task of calculating the expected fluxes is now more demanding, since errors 
in the calculation, or the use of flawed data as input can result in the introduction of 
biases in the estimate of these parameters. 

The study of the zenith angle distributions of the neutrino fluxes is obviously of central 
importance in the determination of the oscillation parameters. The neutrino pathlength L 
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(and the density profile encountered by a neutrino during its path) are strictly correlated 
with the zenith angle 9^ (the correlation is not perfect because the neutrino creation point 
along the straight line defined by its momentum at detection is not exactly known). The 
neutrino flavor transition probabilities depend on the pathlength L (or more in general on 
the neutrino path), and therefore leave their imprints on the zenith angle distributions. 

The azimuth angle distributions are much less interesting. The only significant source 
of non-flatness for these distributions are the effects of the geomagnetic field. In fact, 
neglecting the existence of the magnetic field, the geometry of the neutrino source volume 
(a sherical shell of air) and of the volume where neutrinos propagate, with very good 
approximation have cylindrical symmetry for rotations around the vertical axisQ. Neutrino 
oscillations (or in fact any other forms of transition or disappearance proposed as an 
explanation for the atmospheric v data) do not disturb this symmetry^. The fact that the 
azimuth angle distributions are independent from the new physics that is investigated, 
has however a positive consequence, their study can provide an important cross check for 
the theoretical predictions and the experiments. The calculations must be able to predict 
the shape of these distributions, without the need of additional (unknown) physics beyond 
the standard model, and experiments must be able to measure angular effects that are 
unambiguously predicted. 



1.1 The hint of a discrepancy between data and prediction 

The fluxes of (positively charged) primary cosmic rays that reach the Earth's atmosphere 
contain more particles traveling from west toward east than in the opposite direction. The 
discovery of this effect in 1933 P, |^ allowed to determine that the dominant component 
of the cosmic rays is positively charged. This "primary asymmetry" is also reflected in 
the intensity of the fluxes of secondary particles that are generated in the showers of the 
primary particles and have a direction correlated with the primary one. 

The Super-Kamiokande (SK) detector has measured the azimuth angle distributions 
of atmospheric neutrinos, finding as expected that the distributions are not flat, and that 
there in an excess of particles traveling toward east. In the SK analysis the east-west 
asymmetry is defined as: 

{Ne - Nw) , . 

where A^^; {Ny/) is the number of events with the detected charged lepton traveling toward 
east (west). After selecting events with the charged lepton (£ = e,/i) in the momentum 
interval ^^(GeV) = [0.4,3] and the zenith angle region |cos6'£| < 0.5, the measured 
asymmetries are: 

Af^ = 0.21 ± 0.04, Af^ = 0.08 ± 0.04 (2) 
This has to be compared with predictions from Honda, Kasahara, Kajita and Midorikawa 

^ Other sources of asymmetry can in principle be the presence of mountains (especially above the detector) and the 
existence of different profiles of the atmosphere in different geographycal locations. The first effect is easily calculable. Both 
effects are expected to be negligibly small. 

^ If an east— weast asymmetry is generated by some other source (for example geomagnetic effects), neutrino oscillations 
can in principle affect the size of the asymmetry, since oscillations can distort the u energy spectrum. This effect is of 
second order, and a small correction to the calculated cast- west asymmetry. 
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A™ = 0.13 A™ = 0.11 (3) 

^Bartol _ Q_^7 ^Bartol _ g^^g 

Considering the statistical errors and the spread in the theoretical prediction, the SK 
collaboration concludes that both measured asymmetries are compatible with the expec- 
tations. However, if we consider the ratio Ae/A^, we can see that there is a hint (at the 
level of 2-2.5 standard deviations) of a discrepancy between the calculations that predict 
asymmetries that are approximately equal for e-like and /x-like events (with only a small 
excess for e-like events) and the data where Ag is events is more than twice larger that 
the asummetry for /x-like events. We also note that there is an indication that the mea- 
sured asymmetry for e-like events is larger than the prediction while on the contrary the 
prediction for /i-like events is smaller. 

A discrepancy between data and calculation for the azimuthal distributions can have 
a non negligible importance a for the interpretation of the data. The study of the east- 
west aymmetry probes the atmospheric neutrino fluxes in the angular region around the 
horizontal plane that is actually the most important one for the determination of Am^. 
The fact that the results for e-like and /i-like events differ from the expectations, and 
in opposite directions, in a situation where u oscillations cannot play a role, is clearly 
important for the studies of neutrino oscillations, where the comparison of the events 
rates for the two lepton types is a cornerstone of the analysis. 

In this work we predict that with increasing statistics the Super-Kamiokande mea- 
surement of the east-west asymmetry will clearly show to be in disagreements with the 
predictions of ^ and and indeed with any other prediction based on a calcula- 



tion that does not take into account the effects the bending in the geomagnetic field of 
secondary charged particles in the cosmic ray showers (even if the calculation is three di- 
mensional ||lll)- The existing calculations predict an east-west asymmetry that is simply 



the reflection of the asymmetry of the primary cosmic ray. In this case the asymmetries 
for e-like and /x-like events (and indeed for all for neutrino types Ue, T^e, J^^t, J^/^)- As we 
will discuss in the following the effect of the magnetic field in the development of the 
shower in the atmosphere results in different asymmetries for the different neutrino types, 
in good agreement with the Super-Kamiokande data. 

This work is organized as follows. In the next section we rapidly summarize the effects 
of the geomagnetic field on the primary cosmic rays. Section 3 contains a discussion of 
the effects of the geomagnetic field on the shower development. Section 4 discusses the 
effects of the bending of the muon trajectories in in the geomagnetic field. This is the key 
mechanism that enhances or suppresses the asymmetry for the different neutrino flavors. 
In section 5 we show the results of a complete calculation that includes the bending in 
the geomagnetic fields of all secondary charged particles. Section 6 gives a summary and 
some conclusions. An appendix contains a brief discussion on the difficulty of performing a 
detailed and accurate "three-dimensional" calculation of the atmospheric neutrino fluxes 
and on the possible strategies that are under study. 
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2 The East— West effect for the primary cosmic rays 

The flux of primary cosmic rays protons and nuclei that arrive in the vicinity of the Earth 
surface exhibits the well known east-west effect, that is there are more particles traveling 
from west toward east than in the opposite direction. This is due to the effects of the 
geomagnetic field that forbids the lowest rigidity particles from reaching the Earth. The 
effect is strongest for west-ward going positively charged particles. 

As an illustration let us approximate the geomagnetic field as a dipole, and study the 
trajectories of charged particles in the equatorial plane. The magnetic field is orthogonal 
to the plane, pointing "up" toward the magnetic pole in the northern hemisphere. Near 
the surface of the Earth the field has an approximately constant value B ~ 0.31 Gauss. A 
positively (negatively) charged particle with unit charge and momentum p ~ 59 GeV can 
have a a clockwise (counterclockwise) circular trajectory that grazes the Earth along the 
equator going from east to west (from west to east)0. It is clear that an observer at the 
equator looking near the horizon toward magnetic east (west) can only see unit charge 
primary particles with a momentum larger than 59 GeV, all lower momentum particles 
have "forbidden trajectories". If the equatorial observer turns around by 180 degrees, 
he/she can see primary particles of much lower momenta. This is qualitatively the origin 
of the east-west effect. 

The problem of calculating the geomagnetic effects on the primary cosmic rays (before 
they reach the atmosphere) can be reduced to the problem of calculating "allowed" and 
"forbidden trajectories". The definition of allowed and forbidden trajectories is the fol- 
lowing. Let us consider a cosmic ray of charge q, momentum p near the Earth surface in 
position X, and let us study the past trajectory of the particle. This study can have three 
results: 

(a) the trajectory originates from the Earth's surface; 

(b) the trajectory remains confined in the volume i?® < r < oo without ever reaching 
"infinity" ; 

(c) the particle in the past was at very large distances from the Earth. 

Trajectories belonging to the classes (a) and (b) are considered as "forbidden" , because 
no primary cosmic ray particle can reach the Earth from a large distance traveling along 
one of these trajectories. All other trajectories are allowed. 

The algorithm (used in all calculations of the atmospheric neutrino fluxes) to obtain 
the primary cosmic ray fluxes as a function of position and angle is based on two steps: 

1. It is assumed that the primary cosmic rays at a distance of one astronomical unit 
from the sun are isotropic, unless they are disturbed by the presence of the Earth. 

2. The flux arriving at an imaginary surface in the vicinity of the Earth (the "top of 
the atmosphere") is equal to the isotropic primary flux after the subtraction of all 
"forbidden trajectories" . 

^The trajectory is unstable even in the case of an exactly dipolar field 
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The theoretical basis for the "subtraction agorithm" is sohdly motivated |jT2[ , and is based 
on two fundamental assumptions, that the cosmic ray spectrum "in the absence of the 
Earth" is isotropic, and that the field around the Earth is well described by a static purely 
magnetic field. The Liouville theorem states that the density of points in phase space 
volume is constant. The flux x) is in fact proportional to the phase space density for 
relativistic particles, and since the momentum of a charged particle in a magnetic field 
is constant, then the differential flux (f){p,x{t)) along a particle trajectory is constant. It 
follows that if the flux is isotropic at large distances from the Earth, then the flux in a 
small cone around any trajectory is constant and independent from the position. If a 
particle can reach the Earth, then the momentum spectrum, and the angular distribution 
around it are not deformed. 

For the explicit calculation of allowed and forbidden trajectories, the most accurate 
method is the so called "back-tracking method" that is a straightforward direct applica- 
tion of the definition outlined above. Given a detailed map of the magnetic field around 
the Earth (for example see |]13|) it is a straightforward exercise to compute numerically 
the past trajectories of a cosmic ray from a point "just above" the atmosphere, and test 
if it corresponds to an allowed or forbidden trajectory. 

It is well known that the problem of calculating allowed and forbidden trajectories can 
be solved analytically for the special case of a volume that is entirely filled with an exactly 
dipolar magnetic field (there are no class (a) trajectories in this model). In this case the 
geomagnetic effects result in a sharp rigidity cutoff"Q. For a given position x and a given 
direction -h the trajectories of all positively (negatively) particles with rigidity R > Rg 
{R < Rg) are allowed and the trajectories of all particles with R < Rg {R > Rg) are 
forbidden because they remain confined to a finite distance from the dipole center. The 
quantity Rg {x, h) is the Stormer rigidity cutoff : 



Rtir.XM.O.^) = (^) (5) 

\2r'^ / [ [1 + (1 — cos'^ Am smfc'smv^j^/^J^ J 

where we have made use of the cylindrical symmetry of the problem, M is the magnetic 
dipole moment of field, r is the distance from the dipole. Am the magnetic latitude, 
6 the zenith angle and if an azimuth angle measured clockwise from magnetic north. 
The Stormer cutoff has been of considerable historical importance. It is not sufficiently 
accurate for a modern detailed numerical study, but it remains a good approximation 
that contains all qualitative features of an "exact" numerical calculation, and remains a 
very valuable tool to gain physical insight and qualitative understanding. We can use it to 
ilustrate the three most important qualitative points that we will need for the discussion 
in this work. 

1. The geomagnetic effects depend on the detector position. They are strongest near 
the equator and weakest at the magnetic poles. In equation (|^) we can see that the 
value of the cutoff grows monotonically from a vanishing value at the magnetic pole 
(Am = ±90°) to a maximum value at the magnetic equator (Am = 0°). 



*In a numerical calculation for a non exactly dipolar field, or even in the case of a dipole, if the trajectories that have 
a segment inside the Earth are also considered as forbidden, there is not anymore a sharp cutoff, but in a narrow interval 
allowed and forbidden bands of rigidity alternate. 
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2. The set of allowed rigidities is largest for particles traveling toward magnetic east. 
In equation (|^) for a fixed detector position (fixed Aj\/) and a fixed zenith angle 
9, the cutoff is minimum (maximum) for azimuth angle ip = 270° {(f = 90°) that 
corresponds to a particle traveling toward east (west). 

3. For a fixed detector position and a fixed azimuth angle the set of allowed rigidities 
is largest (smallest) for vertical (horizontal) particles. In equation (|^) for Aj\/ and (f 
fixed, the cutoff grows monotonically with 6. 

The east-west aymmetry of the primary cosmic ray radiation is of course reflected 
in asymmetries for the fluxes of secondary particles generated in their showers in the 
atmosphere. This effect results in asymmetries of approximately the same size for all 
four neutrino types. This is a consequence of the fact that a primary particle of energy 
Eq produces yields of z/g, Ve, and with approximately the same energy spectrum 
(and relative normalizations in the ratios 1:1:2:2). The existing calculations of 
the neutrino fluxes only include the geomagnetic effects on the the primary flux as a 
source of an azimuthal asymmetry and as a consequence the predicted east -west neutrino 
asymmetries are all of approximately the same size. 

3 Effects of the geomagnetic field on the shower development 

3.1 Bending of the trajectories of charged particles 

The trajectories of all charged particles in cosmic rays showers are curved because of the 
presence of the geomagnetic field. It will be particulary important in this work to consider 
the effect of the bending on the zenith angle of the particles, that is to study if the particle 
trajectories are bent "upward" or "downward". We can define a (position dependent) 
system of coordinates with the z axis pointing up, the x axis pointing toward magnetic 
north and the y axis completing a right-handed system pointing toward the magnetic 
west. In this system the magnetic field has by definition components B = (iShorl, 0, -Bvert) 
with invert positive (negative) in the southern (northern) hemisphere. The equation of 
motion of a charged particle in the geomagnetic field for the vertical component is: 

= —qPy l-Bhorl = +5'/^cast |-Bhor|, (6) 

where q is the electric charge and (3 the velocity of the particle. Positively charged 
particles {q > 0) traveling east {(3y < 0) have dpz/dt > 0, and are bent up. Negatively 
charged particles with the same initial direction are are bent down. The reverse happens 
for particle traveling in the opposite (west-ward) direction {Py > 0): negatively charged 
particle are bent up while positively charged ones are bent down. 

An example of the bending of charged particles is shown in fig. |l|. The figure is in true 
scale and represents a projection of the space above the Super-Kamiokande detector. The 
center of the SK detector (latitude 36.42° N, longitude 137.31° E, and altitude 371.8 meters 
above sea level ||^) corresponds in the figuere to the point with coordinates (0,0.3718). 
The thick line represents the "sea level" surface of the Earth in the vicinity of SK detector. 
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Figure 1: Deviation of a muon in the geomagnetic field above the SK detector (see text). 



The z axis corresponds to the vertical axis that passes through the SK center, and the 
horizontal axis points to geographical west. The thin solid lines describe the projections in 
the plane of the figure of the trajectories of three charged particles. The particles are part 
of a cosmic ray shower and are "connected" by a production chain p — 7r+ ^ The 
trajectories are calculated integrating the equations of motion for charged particles in the 
International Geomagnetic Reference Field (IGRF [0) for the year 2000 and correspond 
to a 10 GeV proton, a 7 GeV and a 5 GeV /x"*". The pion production and decay 
points are marked by two diamonds. The muon emits a neutrino (the trajectory and its 
continuation beyond is indicated by a dotted line) that intersects the SK detector. For 
illustration both the pion and the muon decay after exactly three proper lifetimes. Both 
the pion and the muon are emitted exactly collinearly with respect to the primary (or 
parent) particle, so the angle between the neutrino and the primary directions is entirely 
due to the effects of the geomagnetic field. Since all three particle are positive and travel 
toward the east direction they are bent upward. 

In the following we will study the effects of the bending of the shower particles in the 
magnetic field on the fiuxes of atmospheric neutrinos. These effects have been ignored 
in all previous calculations. In our numerical work we have included the bending of all 
charged particles in the shower. In the qualitative discussion that follows we will discuss 
separately the effect of the binding on the primary particles, and on the muons produced 
in the showers. These are in fact the most important sources of the effects we have found. 
Charged pions have a a lifetime that is two orders of magnitude shorter than (and a 
mass similar to m^), therefore their bending in the magnetic field is much less important 
than for muons. Most of the neutrinos are produced in the chain decay of mesons produced 
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in the first interaction of a primary. Also for those produced deeper in the shower most of 
the interesting effects is the result of the bending of the primary particles and the muons. 

Note that we can expect that the effects of the field will be especially important for 
the neutrinos that come from muon decay, and significantly smaller for the muon (anti)- 
neutrinos produced directly in meson decay. 

3.2 The bending of the primary particles 

It is incorrect to think that all effects of the magnetic field on the primary particles 
are taken into account by the study of allowed and forbidden trajectories discussed in 
section 2. The bending of the primary particles is a factor in determining the average 
altitude of the primary interactions, and the zenith angle distribution of the primaries 
at the interaction point. These two effects are not taken into account in the existing 
calculations. Note that the zenith angle (that is the polar angle with respect the local 
vertical) of a particle can be calculated as: 



and is not a constant even for particles traveling along a straight line. The value of the 
zenith angle along a particle trajectory will of course be determined by the bending in 
the magnetic field. 

As an illustration let us study the trajectories of two protons both on allowed tra- 
jectories that cross an imaginary spherical surface above the Earth with the same zenith 
angle but different azimuth angles, with one particle traveling toward east, and second one 
toward west. If the two trajectories are approximated as straight lines the distributions 
of altitude of the interaction point for the two particles (determined by the interaction 
length Aint, the inclination of the trajectory and the altitude profile of the air density) will 
of course be identical. It is however easy to see that since the east-going particle is bent 
"up", and the west-going one is bent "down", the first particle will on average interact 
interact higher, and more "horizontally" : 



Both the inequalities (^ and (^ result in an enhancement of the east-west effect for 
the neutrinos. The average number of neutrinos produced in a shower is larger when the 
primary particle interacts higher because the muons have a longer pathlength for decay, 
and the shower develops in a medium of lower density where the decays of mesons is 
enhanced, and muons lose less energy before decay. For the same reasons more inclined 
showers produce more neutrinos. 

In summary we can expect that the bending of the primary particle inside the atmo- 
sphere result in an enhancement of the neutrino east-west asymmetry with respect to 
the estimates that only consider the bending in vacuum. This enhacement will be ap- 
proximately equal for all four neutrino types, it is therefore an interesting effect, but it 
is cannot be the explanation for the (hint of) discrepancy between data and prediction 
represented by equations @, (H) and (||). 



cos 6* 



— * — * 

p ■ r 



(7) 





{hmt)E > {hint)w, 



(8) 
(9) 



8 



3.3 The angle between the neutrino and the primary particle 

For the effects that we want to consider it is essential to consider the difference in direction 
between a neutrino and the primary particle that generated the shower where it was 
produced. We introduce the notation: = Qq Q) Qq^ where is the direction of the 
detected neutrino, Qq the direction of the primary particle, and f2oi/ the angle between 
the neutrino and the primary. The angle Qq^, is the result of the combination of several 
processes. For the neutrinos that are produced directly in a meson decay we can write: 

^^0. = ^0. © ^nu (10) 

where JIott the direction of emission of the parent pion in the primary interaction, and 
fl-rtu the direction of emission of the neutrino in the pion decay (this qualitative discussion 
obviously applies also for kaons). Similarly for neutrinos emitted in the chain decay 
TT ^ we have: 

^Ou = ^On © ^nf, © © (1 1) 

where Q^^^ {^fj,u) is the emission direction of the muon (neutrino) in the decay of the pion 
(muon), and fi^^ is the deviation of the muon in the geomagnetic field (the deviations 
in the magnetic field of the other charged particles are typically two orders of magnitude 
smaller; they are included in our numerical work). 

The average deviation (considered as a space angle) aj for the process j can be easily 
estimated in first order: 

{Pl.)n Mpi.). 5.2° 

ttOTT ~ TTT-TTV' y^^) 

Ptt Pu pA^^^) 

where we have used the approximation ~ Ap^, (since a charged pion energy is approxi- 
mately shared equally between three neutrinos and a e^) and have assumed a transverse 
momentum {p±) ~ 350 MeV; 

p* 1.7° 
p^ p^GeV) 



where p* is the center of mass momentum of the final state particles in the decay tt 

p* 3p* 0.6° 
P,i Pu PuiGeV) 



this is approximately 3 times smaller that (13) because the muons have on average three 



times the energy of their neutrino decay products; 

mJ3 2.0° 

and for the deviation of the muons in the magnetic field: 

(16) 
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where is the muon decay length, i?^ the muon gyro-radius in the geomagnetic field, 
and B the value of the field. In this estimate we have neglected the muon energy loss and 
we have also assumed that the muons decay before hitting the ground. 

There are two fundamental differences between the effect of the magnetic bending 
and the other contributions to Qq^. The first difference is that all deviations, except the 
magnetic bending one scale oc reflecting the forward boost in the interaction or decay 
of relativistic particles. The deviation of the muons in the geomagnetic fleld is (in first 
order) independent from the muon momentum, because both the magnetic rigidity and 
the decay pathlength increase proportionally to p^. Higher momentum muon bend less 
in the field, but live longer and the field can act for a longer time. 

The second and most important difference is that all deviations (except the bending 
one) are azimuthally symmetric, therefore if we consider the deviation in any plane, the 
average value vanishes. The deviation in the magnetic field is of course not an "average 
value" , all muon of the same initial momentum have exactly the same deviation, and the 
deviation is not azimuthally symmmetric but happens in a well defined plane. If we for 
example consider the average zenith angle of the neutrinos produced in the showers gener- 
ated by primaries of fixed direction Qq = (cos6'o, ^p) we find that the average contribution 
of all sources vanishes, except for the effect of the bending in the field. 

{0u) = 0o + {e,B) (17) 

This effect plays a crucial role in determining the neutrino east-west asymmetry. 

4 Magnetic bending of /i^ and the u east— west asymmetries 

The effect of the geomagnetic field on the trajectory of charged particles (equation (|D) 
implies (a scheme of the results is also shown in fig. |^): 

(a) Positively charged particle (g > 0) traveling toward east {(3y < 0) are bent "up" . 

(b) Positively charged particle (g > 0) traveling toward west {(3y > 0) are bent "down". 

(c) Negatively charged particle (g > 0) traveling toward east {Py < 0) are bent "down". 

(d) Negatively charged particle (g > 0) traveling toward west {Py > 0) are bent "up". 

From equation (|l^ we can then deduce that neutrinos produced in the decay of positive 
(negative) muons in showers of east-ward (west-ward) traveling primaries will be on 
average more "horizontal" than the primary particle; conversely the neutrinos produced 
in the decay of positive (negative) muons in showers of west-ward (east-ward) traveling 
primaries will be on average more "vertical" than the primary particle. 

The key remark is now, that this systematic difference between the zenith angle of the 
neutrino and the primary results for the cases (a) and (d) in an enhancement and for the 
cases (b) and (c) in a suppression of the neutrino flux. 

Two effects play a role in the enhancement or suppression of the flux. The flrst effect is 
purely geometrical and would be present even for an exactly isotropic primary flux (this 
is not so academic as it may sound, because the primary flux is actually isotropic for all 
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Figure 2: Relation between the directions of of the primary particle and of the detected neutrinos, for 
i^'s produced in muon decay (see text). 



rigidities above 60 GV). In this case the rate of primary interactions with zenith angle 
Qq is oc cos6'o, and there is an excess of vertical showers over horizontal ones. Therefore 
if the zenith angle of the primary particles is smaller (larger) than the neutrino one we 
can expect a larger (smaller) event rate. This is illustated schematically in the diagram 
of fig. 0. The neutrinos detected by an observer (indicated by a diamond) located on the 
surface of the Earth (shown as a thick line) looking in the solid angle (ifiwest (indicated by 
a dashed cone) are produced in a "patch" (outlined with a thick line) of atmosphere (that 
in the diagram is represented by a thin spherical shell) subtended by the solid angle d£l. 
It can be seen from the figure that this patch of atmosphere "looks" larger (it has a larger 
projected area) for the primary particles if the production chain is p — vr"*" — >^ /z"*" — >■ 
(the neutrinos can be v^, and than in the case when the neutrinos are produced directly 
in meson decay. This is a simple consequence of the fact that the projected area of the 
element of atmosphere subtended by the solid angle d^l is proportional to cos6'o. For a 
fixed neutrino zenith angle, the average value of the primary zenith angle is not a constant 
but depends on the production mechanism. In the case p — ^ 7r+ — /i"*" — > z/ the relation 
between the zenith angle of the primary and of the neutrinos is {Qy) < 6q, in the case 
of neutrinos directly from meson decay: {O^,) ~ ^o- Similarly the patch of atmosphere 
"looks" smaller for the primaries when the production process is p ^ tt^ ^ fi~ ^ u 
(the neutrinos are Ve and z/^) since in this case the relation between the zenith angle 
of the neutrinos and the primary is: {6^) > Oq. The bending in different directions of 
positive and negative muons results in one case in an enhancement, and in the other in 
a suppression of the neutrino fluxes. This discussion can be repeated for a solid angle 
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(if^east in the opposite hemisphere, with the crucial difference that the production chains 
that are enhanced and suppressed are interchanged. 

If now we consider the fact that the primary cosmic rays that reach the Earth's at- 
mosphere are not isotropic we find actually that the enhancement and suppression are 
stronger that what is obtained with the geometrical argument outlined above. This can 
be uderstood noting that the set of set of forbidden trajectories becomes larger when the 
zenith angle increases. This can be easily seen from the Stormer formula (^) where for 
fixed \m and if the rigidity cutoff grows monotonically with the zenith angle 9. 

In summary: the effect of the magnetic bending of the muons results in an enhancement 
for the flux of the neutrinos produced in the decay of /x"*" (that is Ve and V^) when the 
particles travel toward east, and a suppression when the particles travel toward west. The 
net effect is a positive contribution to the east-west asymmmetry. The effect of magnetic 
bending is obviously opposite for the flux of neutrinos produced in the decay of ^~ {v^ 
and v^). The flux is suppressed when the when the particles travel toward east, and 
enhanced when the particles travel toward west. The net effect is a negative contribution 
to the east-west asymmmetry. 

The effect of the magnetic bending of the muons has to be combined with the asym- 
metry generated by the geomagnetic effects on the primary particles, that to a good 
approximation result in an equal contribution to the asymmetry of all neutrino types. 
Moreover, for the neutrinos produced directly in meson decay, the effects of the geo- 
magnetic field during the development the shower are negligible, and for them the only 
significant source of asymmetry are the geomagnetic effects on the primary trajectory. 

4.1 Qualitative predictions 

We can finally put all results together and give some simple, non-trivial and testable 
predictions. 

1. With increasing statistics, the Super-Kamiokande detector will detect east-west 
asymmetries that are not in agreement with predictions based on one-dimensional 
calculations. 

2. The asymmetries for the four neutrino types are in the relation: 

A^^ < A,^ < AID < ^ ^^g) 

where A^P is the asymmetry predicted in a one-dimensional approximation, that is 
approximately equal for all neutrino types. 

3. The neutrino type with the largest asymmetry is the Ug. The majority of these 
neutrinos are produced in the chain decay p tt^ ^ fi~^ u^, and for all these 
neutrinos the contributions of the geomagnetic effects on the primary particle and 
on the muon add to each other. 

4. The east-west asymmetry for is smaller that for A^^ but larger than the 1-D 
prediction. This is a consequence of the fact that approximately one half of the F^'s 
are produced in /x"*" and have a large asymmetry (approximately equal to A^J and for 
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these neutrinos the geomagnetic effects on the primary particles and on the muons 
add to each other. The other half of the F^'s is produced in pion decay, and will 
have an asymmetry close to what is predicted in the 1-D approximation. 

5. Developing a similar argument we can conclude that the asymmetry for z/^'s will be 
smaller that the 1-D prediction, because for approximately half of these neutrinos 
(those produced in the chain decay p — vr" — — > z/^) the effects of the bending 
of the muons gives a negative contribution to the asymmetry. 

6. The asymmetry for Vf, will be the smallest one since most of the particles are produced 
in the decay of negative muons, and for nearly all Fe's the effect of the bending of 
the muons must be subtracted from the asymmetry generated by geomagnetic effects 
on the primary particles. Numerically we actually will find a deficit of east-going 
Fg's, that is an asymmetry that has changed sign. The effect of the muon bending 
overcompensates for the effect of the primary flux. 

7. Even if the detector is not capable to distinguish neutrinos from anti-neutrinos, the 
effect of the magnetic bending of the muons is still measurable. For e-like events the 
net effect is an enhancement of the asymmetry with respect to the 1-D expectation: 

A, > Af" (19) 

This can be simply deduced as a consequence of three steps: (i) the effect of muon 
magnetic bending on the asymmetry is positive for u^s and negative for Fg's; (ii) 
in the neutrino flux the Ue/Vg ratio is approximately unity (more accurately ~ 1.2, 
reflecting the tt^/tt^ ratio in the final states of proton interactions); (iii) the cross 
section for neutrinos is larger than for anti-neutrinos. 

8. With a very similar argument, we can estimate that the asymmetry of /i-like events 
must be smaller than the 1-D calculation: 

A, < (20) 

The inequality is opposite with respect to the e-like case, since the effect of the muon 
magnetic bending on the asymmetry is negative for neutrinos (with the larger cross 
section) and positive for anti-neutrinos. 



5 Numerical results from a full 3— D calculation 

To go beyond this qualitative discussion we have performed a full three-dimensional mon- 
tecarlo calculation of the atmospheric neutrino fluxes, including the bending in the ge- 
omagnetic field of all charged particles in the cosmic ray showers. This montecarlo cal- 
culation has been already described in [^], where it was used to illustrate effects of a 
three-dimensional calculation on the predicted zenith angle distibutions of e-like and fi- 
like events. We refer to [O for more details on the calculation method, here we will 



only briefly summarize the main points. We have used a straightforward, direct method 
(see the appendix for a comparison with other approaches). The first step of the cal- 
culation is the generation of isotropic fluxes of cosmic rays (we used the results of |]l 
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for the energy spectrum and mass composition) at "the top of the atmosphere" (at a 
radius R = i?® + 80 Km). The geomagnetic effects on the primary flux are computed 
studying the past trajectories of each generated particle. Primary particles on forbidden 
trajectories are rejected. Primary particles on allowed trajectories are propagated in the 
magnetic field (the International Geomagnetic Reference Field for the year 2000 |T^). 
About 97% of the particles interact in the air (the remaining fraction only grazes the 
Earth and continues its travel). The air is described as having a spherically symmetric 
distribution p(r) obtained from a fit to the US standard atmosphere. The final states of 
the hadronic interactions was generated using the model (and montecarlo algorithms) of 
Hillas |15|. All secondary charged particles are propagated along curved trajectories in 



the magnetic field, and a shower is generated developing a "tree" with a standard tech- 
nique. For the propagation of muons we have included the (crucially important) energy 
loss for ionization in the air (and the trajectory correctly takes into account the variation 
of the muon momentum). Energy loss was neglected for all other particles. Multiple 
scattering was neglected for all particles. When neutrinos are produced their trajectories 
(simple straight lines) are studied. Each neutrino either "misses" the Earth, or crosses its 
surface twice, first as down-going and then as up-going. At each intersection its position 
and direction of the neutrinos is recorded. To compute the azimuth of the neutrinos, for 
each point on the surface of the Earth we have defined a "local" system of coordinates 
as discussed in section 3: the z axis points up, the x axis toward magnetic "north", and 
the y axis toward magnetic "west". In this system the magnetic field has components 
(|i?hor|, 0, -Bvert)- The azimuth ip (following the SK convention) is defined as the direction 
where the particle is going: a particle with (p = 90° (270°) is traveling from east toward 
west (from west toward east). 

The neutrinos generated with these algorithms are distributed over the entire surface 
of the Earth, with important non-uniformities due to the effects of the magnetic field. 
We have collected the neutrinos in five different regions of equal area selected according 
to the magnetic latitude: an equatorial region, two intermediate and two polar regions. 
The results in the two (north and south) polar regions and the two intermediate regions 
are essentially undistinguishable, and here we present only the results for the northern 
regions. 

Our results on the azimuth distributions are collected in 6 figures (from fig. 3 to 
fig. 8). In the figures we plot the azimuth distribution of the event rates, obtained after 
the convolution of the calculated fluxes with the model of the neutrino cross section 



from [0. We have selected neutrinos in the momentum interval: p^{GeV) = [0.5,3]) 
and the zenith angle region cos^i, = [—0.5,0.5]. No detection efficiency or experimental 
smearing has been included. For each of the three regions of the Earth that we have 
considered, we include two figures. The first figure contains (in four separate panels) 
the azimuth distributions of the four neutrino types. The second figure (in two separate 
panels) contains the azimuth distribution of e-like and /x-like events (these distributions 
are obtained simply summing together the results for and Ve (or and u^) . The scale 
of the y axis is absolute for all figures. 

We have also performed the montecarlo calculation three times. 

(i) A first calculation (represented by thick histograms) was performed using the "full 
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3-D" algorithms descrivbed above, including the bending of all secondary charged 
particles in the showers. 



(ii) A second calculation (represented by thin histograms) was performed to reproduce 
the 1-D algorithms. The geomagnetic effects are calculated for the primary cosmic 

ray particles traveling outside the Earth's atmosphere (exactly as in the previous 
case), but all particles travel along straight line trajectories for r < + 80 Km.). 
Also all final state particles are coUinear with the projectile (or parent) particle. 
This is achieved modeling the interactions and particle decays exactly as in the 
previous case, (including therefore transverse momentum), and performing as a last 
step a rotation of all the 3-momenta of the final state particles so that they become 
parallel to the projectile (for interactions) or parent (for decays) particle. 

(iii) A third calculation (represented by dashed histograms) was performed neglecting 
the geomagnetic effects on the primary flux (therefore considering exactly isotropic 
primary fluxes) and using the 1-D algorithms outlined in the previous point. With 
these approximations the azimuth angle distributions must be flat and independent 
from the detector position. The only non trivial result of the calculation is the 
absolute normalization of the different rates. Note that this normalization must also 
be independent from the geographycal region considered. 

Inspecting the results of flg. 3 to flg. 8, we can make the following remarks. 

1. Calculation (iii) (no geomagnetic effects) results as expected in flat, detector position 
independent azimuth distributions. 

2. Calculation (ii) (ID with geomagnetic effects included for the primary) exhibits as 
expected an east-west asymmetry, that is of approximately the same size for all four 
neutrino types. 

3. The taking into account of the effects of the bending in the magnetic field of secondary 
particles (and also of the primary in the atmosphere) in calculation (iii) is a non 
negligible correction. 

4. Comparing the results of calculation (ii) and calculation (iii) one can see that the 
effects of the bending of secondary particles result in an enhancement of the east-west 
asymmetry for i/g, and 17^, and a suppression of the asymmetry for I7g and f^. 

5. The asymmetry for the combination (i/g + i^e) is also enhanced, while the asymmetry 
for the combination (i/^ + 17^) is suppressed. 

In summary we can observe that the expectations of the qualitative discussion of the 
previous section have been confirmed in a preliminary, but detailed, quantitative calcu- 
lation. The results are summarized in table 1, that gives the results of the asymmetries 
obtained with the calculation (ii) (or 1-D) and the calculation (iii) (or 3-D), in the dif- 
ferent geographyical regions. For the definition of the magnetic latitude we have used as 
dipole axis the one that corresponds to the leading term in the multipole expansion of 
the IGRF field for the year 2000. With this definition the SK detector has a magnetic 
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Table 1: East-west asymmetry for different event types, after averaging for the detector position in 
different regions of the Earth. The events are selected with the cuts Ei,{GeV) — [0.5,3] and cos 9^ = 
[—0.5,0.5]. The two results are for a one-dimensional and a three-dimensional calculation that includes 
the bending of secondary particles in the geomagnetic field. 



Region 


Equatorial 


Intermediate 


Polar 




sin Am = 


= [-0.2,0.2] 


sIuAa/ = 


: [0.2,0.6] 


sin Am = 


= [0.6, 1] 


Event type 


A (ID) 


A (3D) 


A (ID) 


A (3D) 


^(ID) 


A (3D) 




0.194 


0.431 


0.126 


0.335 


0.026 


0.214 




0.189 


-0.057 


0.120 


-0.065 


0.012 


-0.061 




0.180 


0.065 


0.115 


0.028 


0.011 


0.002 




0.173 


0.310 


0.107 


0.240 


0.014 


0.153 


e 


0.192 


0.301 


0.124 


0.224 


0.022 


0.136 


M 


0.178 


0.139 


0.113 


0.091 


0.012 


0.046 



latitude of 27.08° (sin = 0.455) and is inside the "intermediate" region in the northern 
hemisphere. 

Inspecting table 1, we can notice that as expected the east-west asymmetry depends on 
the magnetic latitude of the detector, and is strongest near the magnetic equator, where 
all geomagnetic effects are enhanced. In the 1-D calculation the asymmetries for the 
four neutrino types are approximately equal, with small differences that reflect the fact 
that the neutrino spectra produced by a primary of fixed energy are not identical for the 
four flavors. The full 3-D calculation clearly exhibits the enhancements and suppressions 
predicted qualitatively in section ^. 

6 Summary and conclusions 

In this work we have shown that the azimuth angle distributions of atmospheric neutrinos 
are shaped by the effects of the geomagnetic field on both the trajectories of the primary 
cosmic rays, and their showers in the Earth's atmosphere. 

In the existing (one-dimensional) calculations of the atmospheric neutrino fluxes the 
effects of the magnetic field on the shower development are ignored. In a 1-D framework 
the asymmetries of the v fluxes reflect only the highest suppression of the primary cosmic 
ray flux when the primary particles travel towards west, and the predicted east -west 
asymmetry is approximately equal for all four neutrino types: (z/g, V^, z/^ and V^). 

The situation changes when one includes the bending in the geomagnetic field of sec- 
ondary particles (and in particular of muons) in the cosmic rays showers. Positively and 
negatively charged particles are curved in different directions and the resulting effect on 
the east-west asymmetry is different for different u types. It is an enhancement for the 
neutrinos that are the product of jj,^ decay (z/g, and V^) and a suppression for the products 
of /i~ decay (z/g, and z/^). 

For detectors like Super-Kamiokande that cannot separate z/'s and z7's the detectable 
effect is an enhanced east -west asymmetry for the e-like events, and a suppressed asym- 
metry for the /i-like events, as it is already observed |]^ by Super-Kamiokande. It is 
remarkable that the detctor is able to measure successfully subtle effects as those de- 
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scribed here (and without the need of a prediction). In our view this is a confirmation of 
the high quahty of the experimental data. 

The (partial) failure of the existing calculations in the prediction of the east-west 
asymmetry is also a warning and an indication that the development of more accurate 
and detailed calculations of the atmospheric neutrino fluxes is needed. Several groups 
[|T^] are working in this direction. We note that a correct prediction of the east-west 
effects is possible only with a fully three-dimensional calculation. In our view the success 
of the calculation presented here in reproducing the SK results is also a confirmation of 
the non-trivial geometrical effects connected with 3-D effects and the geometry of the 
neutrino source volume discussed in |Tl| and in more detail in fl^ . 

The understanding that the present calculations of atmospheric neutrinos do not de- 
scribe correctly the east-west asymmetry can have some interesting consequences for the 
determination of the oscillation parameters from the atmospheric neutrino data (in SK 
and also other experiments, in particular Soudan). We note that the effects we have been 
discussing are obviously especially important for horizontal neutrinos, and that the effects 
have different sign for e-like and /i-like events. The shape of the zenith angle distribu- 
tion of the atmospheric neutrino fluxes, in particular close to the horizontal plane, and 
the relation between the muon and electron event rates are (at least in the opinion of 
this author) the most important problems in the calculation of the atmospheric neutrino 
fluxes. As an illustration, in the simplest (and favored) picture of ^ Vr oscillations for 
values of Am? in the range indicated by the data and the range of E^, characteristic of 
the contained events, the value of sin^ 20 can be reliably calculated comparing the "up- 
going" and "down-going rates", with a relatively small systematic error and even for a 
poor (or even incorrect) estimate of Am?. This is possible since the ratio of the rates 
is to a good approximation determined only the "asymptotic forms" of the oscillation 
probability {Py = or sin^2^/2) where the value of Am^ is absent. On the other 
hand the determination of Am^ requires a detailed fit of the zenith angle distributions, 
with the region close to the horizontal playing a crucial role. The (nearly exact) up-down 
symmetry of the z/-fiuxes is of little help in this case. The comparison of the e-like and 
/i-like rates remains as a powerful tool to estimate the suppression (or enhancement) of 
a flavor type, but of course as discussed in this work, for precise quantitative evaluations 
of the parameters one need to study with great care the relations between the fluxes of 
different neutrino types. 

We postpone to a future work a quantitative estimate of the effects of a full three- 
dimensional calculation in the determination of the allowed region in the neutrino oscil- 
lation parameter space. 



Other important sources of uncertainties pO| , pT| in the determination of the allowed 



region for the oscillation parameters are the input primary spectra, where very valuable 
new data became recently available p2|, ESl , and the modeling of hadronic interactions. 



New data on this problem would also be of great value |2^. A detailed description of the 
neutrino cross section is also needed to compute the event rates, and new data would be 
very valuable and benefit also the long-baseline programs. 
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Appendix: 3— D calculations of atmospheric neutrinos 



The first calculations of the atmospheric neutrino fluxes have been performed in a one- 
dimensional approximation, that is assuming that the neutrinos are emitted collinearly 
with the primary particle. This allows an enormous reduction of the size of a montecarlo 
calculation, because only the (formally vanishingly small) fraction of the cosmic rays 
showers that has trajectories that intersect the detector has to be simulated. 

The difficulty of a a full three-dimensional calculation of the neutrino fluxes has been 



discussed in [|lT| and |T3- In a 3-D calculation any shower can produce a neutrino that 
intersect the detector we are considering, and therefore all possible showers have to be 
studied. Since the atmospheric neutrinos are generated quasi-uniformally over the entire 
surface of the Earth with a total area 5.1 x 10^ Km^, only a very small fraction of these 
neutrinos will intersect a dector with an area of order 10~^ Km^. Of course it is possible 
to consider a detector area greatly enlarged to improve the statistical precision of the 
calculation, however this cannot be done arbitrarily because in fact the v flux is not 
exactly uniform, but it does depend on the detector location. This in a nutshell is the 
computational problem of a 3-D calculation. 

There are several methods that have been used or are currently under study to overcome 
this difficulty. 



6.1 Spherically symmetric problem 

A useful first step is to consider a a simpler, sperically symmetric problem, that is obtained 
neglecting all effects of the geomagnetic field both for the determination of allowed and 
forbiddden trajectories, and in the development of the showers. This problem is explicitely 
spherically symmetric: all points on the surface of the Earth are equivalent, and the 
entire surface can be used as the "detector" . This calculation requires therefore the same 
computer power of a 1-D calculation, because essentially all produced neutrinos (all those 
that intersect the Earth's surface) can be collected and analysed to "measure" the i/-fiuxes 
in a montecarlo calculation. 



A calculation along these lines has been performed by G. Battistoni et al |rT|. This 
calculation for the first time demonstrated the existence of non-trivial geomatrical effects 
(see for more discussion) due to the spherical geometry of the neutrino source volume. 



6.2 Shower translation algorithm 

Unfortunately a calculation performed using the approximation of a spherically symmetric 
Earth is not adequate for the analysis of experimental data. The largest effect that is 
missing is the effect of the geomagnetic effect that "forbids" low rigidity cosmic rays from 
reaching the vicity of the Earth. If one neglects the effects of the geomagnetic field on the 
shower development it is possible to "translate" a shower from a position on the Earth to 
an arbitrary one, keeping the zenith angle of the primary particle (and therefore the zenith 
angle of all particles) as constant (this of course corresponds to a rotation of a shower seen 
as "rigid body" around the center of the Earth, plus an additional rotation around the new 
vertical axis). A single shower generated by montecarlo can then be "used" many times 
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with enormous saving of computer time. The effects of the geomagnetic effects in the 
determination of the allowed and forbidden trajectories can be easily included, checking 
that if a given shower position corresponds to an allower of forbidden primary trajectory 
and rejecting forbidden trajectories. A preliminary calculation using this method for the 
three sites of Kamioka, Soudan and Gran Sasso has been recently made available on the 
web by Battistoni et al. PB . 



6.3 Weighted technique 

If the effects of the magnetic fields are also included in the development of the shower, 
the problem loses all symmetry and one has to confront the computational problem of a 
calculation that is intrinsically very inefficient. A solution of this problem that is under 
study ||2^, is to use a weighted technique, generating the cosmic ray showers with a strong 



bias in position and direction, so that the neutrinos produced are more likely to arrive 
in the vicinity of the detector in consideration, and using a weight system to estimate 
correctly the flux. In this method approximately half of the cosmic rays interact in a small 
region above the detector to simulate the down-going flux, while approximately one half 
are generated over the rest of the Earth's surface. For these events the direction of the 
shower is chosen preferentially in a cone with an axis that "points" toward the detector. 
In this way it is possible to collect a sufficiently large statistics of neutrinos in a relatively 
small area around the detector site. 



6.4 Direct approach 

A possible solution is also to neglect any attempt at simplification, or optimization, and 
simply generate all cosmic showers on the Earth with a realistic distributions in posi- 
tion and direction, that only takes into account the geomagnetic effects studying the 
development of the showers in a realistic magnetic field. Such a calculation will produce a 
population of neutrinos that is distributed quasi-uniformally over the entire surface of the 
Earth, with the non-uniformities representing the expected variations in intensity due to 
geomagnetic effects, that produce a larger (smaller) flux in the magnetic polar (equato- 
rial) region. In this direct (or brute force) approach the calculation is not performed for a 
given detector position, because all positions on the Earth are treated "democratically". 



This is the technique used in |jT^ and also in this work. Its only drawback is that the 
results that can be obtained in a reasonable time integrating then neutrino results over a 
large area of the Earth. 
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Figure 3: Azimuth angle distributions of atmospheric neutrino events. The distributions are averaged 
for detectors uniformly distributed on all positions in the Earth magnetic equatorial region (sin Am = 
[—0.2,0.2]). The four panel are for neutrinos of different type. The three histograms are for: fully 3-D 
calculation (thick), 1-D calculation (thin), 1-D without geomagnetic effects (dashed). 
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Region: sin X^^^ =[—0.2, +0.2] 
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Figure 4: Azimuth angle distributions of atmospheric neutrino events. The distributions are averaged 
for detectors uniformly distributed on all positions in the Earth magnetic equatorial region (sin Am = 
[—0.2, 0.2]). The upper (lower) panel is for e-like (/x-like) events. The three histograms are for: fully 3-D 
calculation (thick), 1-D calculation (thin), 1-D without geomagnetic effects (dashed). 
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Region: sin A 
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Figure 5: Azimuth angle distributions of atmospheric neutrino events. The distributions are averaged for 
detectors uniformly distributed on all positions on the Earth with magnetic latitude sin Am = [0.2,0.6]. 
(The SK detector has magnetic latitude sinA|^ = 0.455) The four panel arc for neutrinos of different 
type. The three histograms are for: fully 3-D calculation (thick), 1-D calculation (thin), 1-D without 
geomagnetic effects (dashed). 
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Figure 6: Azimuth angle distributions of atmospheric neutrino events. The distributions are averaged for 
detectors uniformly distributed on all positions on the Earth with magnetic latitude sin Am = [0.2,0.6]. 
(The SK detector has magnetic latitude sinA^^ = 0.455) The upper (lower) panel is for e-like (/i-like) 
events. The three histograms are for: fully 3-D calculation (thick), 1-D calculation (thin), 1-D without 
geomagnetic effects (dashed). 
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Figure 7: Azimuth angle distributions of atmospheric neutrino events. The distributions are averaged 
for detectors uniformly distributed on all positions in the Earth magnetic north polar region (sin Am = 
[0.6,1]). The four panel are for neutrinos of different type. The three histograms are for: fully 3-D 
calculation (thick), 1-D calculation (thin), 1-D without geomagnetic effects (dashed). 
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Figure 8: Azimuth angle distributions of atmospheric neutrino events. The distributions are averaged 
for detectors uniformly distributed on all positions in the Earth magnetic north polar region (sin Am = 
[0.6, 1]). The upper (lower) panel is for e-like (/i-like) events. The three histograms are for: fully 3-D 
calculation (thick), 1-D calculation (thin), 1-D without geomagnetic effects (dashed). 



28 



